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Abstract. Two different types of orthogonality condition models (OCM) are 
equivalently formulated in the Faddeev formalism. One is the OCM which 
uses pairwise orthogonality conditions for the relative motion of clusters, and 
the other is the one which uses the orthogonalizing pseudo-potential method. 
By constructing a redundancy-free T-matrix, one can exactly eliminate the 
redundant components of the total wave function for the harmonic-oscillator 
Pauli-forbidden states, without introducing any limiting procedure. As an 
example, a three-a-particle model interacting via the deep aa potential by 
Buck, Friedrich and Wheatley is investigated. 

1 Introduction 

The Faddeev formalism for three composite particles has always suffered from the 
insufficient treatment of the Pauli principle. For example, in the three-alpha (3a) 
Faddeev study by one of the authors [1] one-term separable aa potentials are used 
to reproduce the damped inner oscillations of the relative wave functions, which 
are the most important effect of the Pauli principle between two a (2a) clusters. 
Due to the insufficient treatment of the Pauli principle among 3a clusters, some 
of the obtained states (02^ and 1" states) were concluded to be spurious since 
they contain a large amount of redundant components. On the other hand, a large 
binding energy of the shell-model like ground state and the excited 0"*" state with 
well-developed cluster structure are simultaneously reproduced, which can never 
be realized by Ali-Bodmer's phenomenological aa potential with the repulsive 
core. The large overbinding of the 3a ground state with £"3^ = — 17 MeV (without 
the Coulomb force) is most easily understood by considering that the damped 
inner oscillations of the 2a relative motion are enhanced in the compact 3a 
system, and the attractive nature in the short-range part of the core- less potential 
overwhelms the large kinetic energies. [2] The same situation is also found in the 
work by Oryu and Kamada [3]. They started directly from the microscopic 2a- 
cluster kernel of the resonating-group method (RGM), and converted it to the 
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many-rank separable potentials which arc suitable for the Faddeev calculations. 
Although their admixture of the redundant components is relatively small, the 
ground-state energy is extremely large with E^a ~ —20 MeV. [4] Since the 
fish-bone optical model proposed by Schmid [5] also gives a large ground-state 
energy, these authors claim that some sort of 3a force is definitely necessary in 
the 3a-particle model to obtain a reasonable agreement with the experimental 
observation. [6] We think that this repulsive 3-cluster force can be partly avoided 
by the complete elimination of the 3a redundant components, which cannot be 
excluded at the 2a potential level. From the microscopic viewpoint based on the 
3a RGM, the model space for the relative motion of the 3a clusters has a well- 
defined notion solely determined from the assumed internal wave functions of 
the a clusters. The standard interpretation of the 3-cluster force may be the one 
which stems from the interaction kernel connected to the full antisymmetrization 
among the three clusters. [7] 

Recently, we have developed 3-cluster Faddeev formalism which employs 2- 
cluster RGM kernel directly. [8, 9] In this formulation, we first write down the 
RGM equation in the form of the Schrodinger-type equation. The resultant in- 
teraction term becomes non-local and energy-dependent. This linear energy de- 
pendence in the interaction term originates from splitting the overlap kernel into 
the direct term and the exchange normalization kernel. The two-cluster RGM 
equation sometimes involves redundant components. In such a case, the com- 
plete off-shell T-matrix is not well-defined in the standard procedure. [10] Our 
strategy is to distinguish between the energy e in the interaction term and the 
starting energy oo involved in the 2-cluster Green function. Assuming the energy 
e involved in the interaction term as a mere parameter, we can define the full 
T-matrix, T(a;, e), through the standard procedure. Although T(a;, e) is singular 
when u! = e, there is no harm in solving the Faddeev equation for the bound 
states since oj is negative and e is usually positive. Our finding is that the modi- 
fied T-matrix, T{u),e), obtained by subtracting this divergent term is the proper 
"RGM T-matrix" , which should be used in the Faddeev equation. The remaining 
energy-dependence in T(a;, e) should be determined self-consistently by calculat- 
ing the expectation value of the 2-cluster Hamiltonian with the resultant total 
wave- function of the Faddeev equation. In Ref. [8], we proved that this formalism 
is completely equivalent to the 3-cluster orthogonality condition model (OCM) 
with pairwise orthogonality conditions (as proposed by Horiuchi [11, 12]), using 
the RGM interaction term as the pairwise interaction. The energy dependence in 
the interaction term should be determined self-consistently even in this 3-cluster 
equation. 

In this paper, we will show that this Faddeev formalism developed for 2- 
cluster RGM kernel is also applicable to the usual 3-cluster OCM. In this case, the 
corresponding "OCM T-matrix" no longer involves the e-dependcncc, although 
the interaction term of the Schrodinger-type OCM equation is still energy depen- 
dent. We will also show that this 3-cluster Faddeev equation is equivalent to the 
Faddeev equation formulated for the "redundancy-free" T-matrix obtained from 
the original OCM potential by applying Kukulin's method of orthogonalizing 
pseudo-potentials [13]. Through these procedures, we can prove the equivalence 
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between Horiuchi-type 3-cluster OCM with the pairwisc orthogonahty conditions 
and the method of orthogonahzing pseudo-potentials. A nice point of the present 
approach is that one no longer needs the process to take the strength parameter, 
A in Eq. (2.30), infinity, which may cause a serious numerical instability if the 
model space to solve the 3-cluster equation is too small in the variational-type 
calculations. We think that this is a great merit of using the T-matrix formal- 
ism. As a typical example, we investigate the 3a system in which the pair a's 
interact via the deep local potential proposed by Buck, Priedrich and Wheatley 
(BFW potential) [14]. According to the general idea of the point-like a-particlc 
models, the redundant components of the 2a system are assumed to be the 
bound states of the 2a Hamiltonian. In this case, the Pauli projection operator 
P needs a careful treatment to select a physical model space. We find that the 
ground-state energy of this system is far below the experimental value. This is a 
different conclusion from that reached in Refs. [15] and [16], which was based on 
the variational calculations in the method of orthogonahzing pseudo-potentials. 

In the next section, the Faddeev formalism for the 3-cluster OCM is developed 
after a brief recapitulation of the previous Faddeev formalism using the 2-cluster 
RGM kernel [8, 9]. An application to the BFW potential with the bound-state 
Pauli-forbidden states is also given, together with a new feature influenced by the 
almost forbidden components of the Faddeev equation. The relationship of the 
T-matrices derived in this particular case is easily generalized for more general 
types of the Pauli-forbidden states composed of the harmonic-oscillator (h.o.) 
wave functions. The third section discusses the numerical examples of the present 
Faddeev formalism for the 3a system interacting via the BFW potential. The 
last section is devoted to a summary. A simple formula for T-matrices is given 
in the Appendix. 

2 Formulation 

2.1 3-cluster Faddeev equation using the 2-cluster RGM kernel 

Wc start from a two-cluster RGM equation for the relative wave function %, 
expressed as 

[u-Ho-V''''^{u;)]x = 0, (2.1) 

where to is the total energy in the center-of-mass (cm.) system, measured from 
the two-cluster threshold, u = E—E™^, Hq is the relative kinetic-energy operator, 
and 

yRGM^^-) = + G + a;if , (2.2) 

is the RGM kernel composed of the direct potential Vd, the sum of the exchange 
kinetic-energy and interaction kernels, G = G^ -\- G^ , and the exchange nor- 
malization kernel K. For simplicity, we assume a single-channel RGM and that 
there exists only one Pauli-forbidden state \u) (normalized as {u\u) = 1), which 
satisfies the eigen- value equation K\u) = 7|u) with the eigen- value 7 = 1. The 
projection operator on the Pauli-allowed space for the relative motion is denoted 
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by yl = 1 — |u)(n|. Using the basic property of the Pauh-forbidden state \u), 

{Ho + Vd + G)\u) = and {u\(Ho + Vd + G) = 0, we find that Eq. (2.1) is 
equivalent to 

A[uj-Ho-V^^^{u)]Ax = (2.3) 

or 

A[uj - Ho - v{uj)] Ax = . (2.4) 

Here, we have defined 

v{s) = AV^'^^ie)A = A{Vb + G)A + eAKA , (2.5) 

by generalizing lu in v(u)) to e. The OCM approximation consists of 

A[iu-Ho-VB]A2P = , (2.6) 

or more favorably changing the direct potential 1^ to a suitable effective local 
potential Vcs- [17, 18] 

The basic procedure to define the "RGM T-matrix" is to separate F^^^^(£) 
into two distinct parts 

V^G^{s) = V{s)+v{£) (2.7) 

with 

V{£) = {e- Ho) - A{£ - Ho) A = e\u) {u\ + AHoA - Ho , (2.8) 

and to assume e in t;(e) as a mere parameter which should be determined by 
the surroundings of the interacting two clusters. Then we can derive a formal 
solution of the T-matrix equation 

T{ij, e) = V^^^ie) + V^^^ie) G^o^\u) T{uj, e) (2.9) 

with G^Q^\uj) = l/{uj-Ho + iO) as follows: [8] 

T{iu, e) = f{iu, e) + {iu - Ho)\u) {u\{u - Ho) , 

UJ — £ 



T{co,e) = T^{co,e)-\l + T^{uj,e)Gl^\co 



u 



{u\Gi-^\iu,£)\u) 
x{u\[l + Gi^^\Lo)mLo,e)j , (2.10) 

where r^(w,£) is defined by 

mLu,e)=vie)+v{e)Gi-^\u)n{u,e) , (2.11) 

and G^\uj,e) is the corresponding full Green function. 

The introduction of the parameter e in Eq. (2.9) is not strange if we consider 
the practical method to solve equations like Eqs. (2.1) and (2.6) with a redundant 
solution \u). In these equations we are actually solving Ax or A^. The original 
Saito's suggestion [17, 18] for solving the OCM equation Eq. (2.6) is to assume 
Alp = ip or {u\ip) = and solve 

cui; = A{Ho + Vb)AiP . (2.12) 
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In fact, the solution of Eq. (2.12) is the redundancy-free solution of Eq. (2.6) when 
a; 7^ 0. Generally speaking, the trivial solution \u) needs not be a zero-energy 
solution. One can also move this exceptional energy to an arbitrary (usually 
positive) value e. This can be achieved by simply adding e\u){u\i/j) term in the 
right-hand side of Eq. (2.12): 

u^ = A{Ho + VB)A'ilj + e\u){u\'ilj) . (2.13) 

One can take the same process in the RGM equation Eq. (2.1). We start from 
Eq. (2.4) and change u in v{u!) to e. This is permissible since the energy de- 
pendence of the v{e) term is usually very weak due to the structure AKA. For 
example, this term vanishes completely for simple systems like the two di-neutron 
system. [8] The energy dependence of the RGM interaction in the allowed space 
is later taken into account by a self-consistency condition. Similarly to Eq. (2.13), 
we set up with the equation 

ux = MHo + v{e)]Ax + e\u){u\x) • (2.14) 
If we use Eq. (2.7), we can easily prove that this equation is nothing but 

[u-Ho-V^^^{e)]x = ^, (2.15) 

which no longer has the trivial solution \u) except for ui = e. 

A motivation to use T(u},e) in Eq. (2.10) for the Faddeev equation comes 
from the complete equivalence between the Faddeev equation and the 3-body 
equation interacting via v{£) in the allowed model space. [8] Namely, for a sys- 
tem composed of three identical spinless particles, we can prove the equivalence 
between 

P[E-Ho- V^'^^isa) - If ^^(£^) - ^^(£^) ] PlF = (2.16) 

and 

V-a = [E) ri^) [E, e^) (V/3 + V7) > (2-17) 
where a common self-consistency condition 

e« = {P^\ K + V^^'^ie^) \P^)/{P^\P^), (2.18) 

is imposed. In Eq. (2.16), Hq is the three-body kinetic energy operator in the cm. 
system, V^'^^{sa) represents the RGM kernel Eq. (2.2) for the a pair, and P is 
the projection operator onto the [3]-symmetric Pauli-allowed space, as defined in 
Refs. [8] and [11]. In the Faddeev^quation Eq. (2.17), Ta {E,ea) is essentially 
the non-singular RGM T-matrix T{u!,£) defined through Eq. (2.10): 

fi^HE,ea) = f^{E-K,e^) , (2.19) 

where ha is the relative kinetic-energy operator between the a-pair and the third 
particle. A nice point of the Faddeev equation Eq. (2.17) is that the total wave 
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function ^, constructed from the three Faddeev components, ipa, ipp and ■^-y is 
automatically orthogonal to the Pauli-forbidden state in each pair: 

{UaW) = {Ua\P'^) = {Ua\^a + # + V7) = ■ (2-20) 

This is because of the orthogonality property 

ti) = , (2.21) 
which is derived from the formal solution in Eq. (2.10). 



2.2 Application to the 3-cluster OCM 

The above discussion on the 3-cluster systems interacting via pairwise 2-cluster 
RGM kernels can be straightforwardly extended to the ordinary 3-cluster OCM 
interacting via simple energy-independent local potentials. One only needs to 
modify 

v{e) — ^v = AVbA . (2.22) 

From Eqs. (2.7) and (2.8), the interaction term for the 2-cluster OCM equation, 
yOCM(-£^^ turns out to be 

V^G^{£) -^V^^^ie) = V{e) + v = {e-Ho)-A{£-Ho-VY))A 

= e\u){u\+A{Ho + VB)A-Ho. (2.23) 

The full T-matrix of V^^^{e) is defined through 

T{uj, £) = V^^^\e) + y°^^(£)Gf,+) (w)T(a;, e) . (2.24) 

The formal expression of T{uj,e) is very similar to Eq. (2.10), but this time 
the "OCM T-matrix", T{uj), docs not involve the e-dependence. Namely, the 
e-dependence appears only through the last term in T(uj,e). Now we can write 
down two equivalent equations corresponding to Eqs. (2.16) and (2.17): 

P[E-Ho-V^ -V^ - F° ]P^ = (2.25) 

and 

V'a = (E) fi^) [E) {^f, + V^) . (2.26) 

This time, we do not need the self-consistency condition Eq. (2.18), and Eq. (2.19) 
becomes 

f^'\E)=f^{E-ha) . (2.27) 

Since the orthogonality condition in Eq. (2.21) is still valid for T{uj), we can 
prove that the solution of Eq. (2.26) satisfies the orthogonality of the total wave 
function, Eq. (2.20). 
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If the Pauli- forbidden state \u) is a real bound state of Vd, our expression for 
yOCM^g) in Eq. (2.23) is further simphfied. We assume that \u) = \ub) is the 
bound-state wave function, satisfying 

(£B - i^o - Vb) = , (2.28) 

with £b (< 0) being the bound-state energy. Here, Hq is the 2-cluster kinetic- 
energy operator, and only one bound state is assumed to exist. Then, one can 
easily show that V*^'^^'^{e) is reduced to 

fOCM(^^^^^^^^_^^)|^^^^^^| ^ (2.29) 
which is nothing but the Kukulin's pseudo-potential 

VP"" = Vb + X\ub){ub\ (2.30) 

with X = e — £b- By using the general formula given in Appendix, we can find 
that the T-matrix defined through Eq. (2.24) is given by 

T{lo, e) = rD(a;) + [ 1 + rD(a;)Go(a;) ] \ub)— 

- {ub\G\)'{uj)\ub) 

x{ub\[1 + Go{u)Tb{u)] , (2.31) 



where 
and 



TBiio) = VB + VBG'^o^\io)TB{uj) , (2.32) 



G^^\u;) = [co-Ho-Vu + iO]-' = Goiu;) + Goiu;)TBiio)Goiu;) . (2.33) 
On the other hand, {u — Hq — Vo)\ub) = {oo — £b)\ub) yields 

g'^^\u:)\ub) = \ub) , (2.34) 

UJ — £b 

if o; 7^ £b- By using Eq. (2.34), the second term of Eq. (2.31) in the right-hand 
side is greatly simplified through the relationship like 

[1 + rD(a;)Go(a;)] l^xs) = {u - Ho)G^+\u)\ub) = {u - Ho)\ub)^—.{2.35) 

u — Eb 



We find 



T{u;,s) = Toiio) -{to- Ho)\ub)^—{ub\{uj - Hq) 

UJ - £B 

+{u - Hq)\ub)-^{ub\{uj - Ho) . (2.36) 

(J — £ 



If we define T{u;) through 



T{u;,£)=T{u) + {u;-Ho)\ub)^—{ub\{uj-Ho) , (2.37) 

LO — £ 
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we obtain 

f (a;) = rD(a;) - (a; - Ho)\uB)^—{uB\{io - Ho) , (2.38) 

U — £b 

or 

f(w) = lim T(w,£) . (2.39) 

Adding the separable term to Vo in Eq. (2.29) removes the bound-state pole 
of the T-matrix Ty){uj) and moves it to the positive energy e in Eq. (2.37). In 
T(ijj) in Eq. (2.38), this positive energy pole is even removed. In order to see 
this, we use the spectral decomposition of G^j^\lj): 

G^j^\u) = [uj-Ho-Vb + ^0]-^ 

Using this in rD(a;) = Vd + Vd g'^\uj) Vd and in Eq. (2.38), we find 
f{L0) = AVbA - \ub) {ub I (w - Ho)\ub) {ub I 

which has no singularities if a; < 0. 

Let us consider a three-body system composed of three identical spinless 
particles with the mass M. One of such examples is the 3a system interacting 
via the BFW potential, discussed in the next section. The pairwise interaction is 
assumed to be the Kukulin's pseudo-potential with A = e — £b in Eq. (2.30). 
The three-body Schrodinger equation reads 

(^E-Ho-VP' -V^' -V:P'^Wr = . (2.42) 

This equation is transformed to the Faddeev equation by the standard procedure: 

ipa = GoTa{tljp + ^P^) , (2.43) 

where Go = G^o^\e) with E < 0, Ta = Ta{E - h^), and etc. arc the three 
Faddeev components yielding the total wave function iZ^ = ■i/^a + V'/J + V't- Using 
Eq. (2.37), we find 

V'a = GoIt^ + {E-Ho)\u^) 3 J {u^\{E-Ho)]{i^p + i^,) 

= GoTc (V/3 + V7) + kg) p._3J {Ua\E - Ho\^p + ■ (2.44) 

Here, q is the momentum Jacobi coordinate between a pair and the third particle. 
If we multiply Eq. (2.44) by {u^\ from the left-hand side, we obtain 
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or 

{Uam = ^ 3^2^ o - ^OIV'/J + V't) ■ (2-46) 

If we take the limit A — > oo or e — > oo, we obtain 

V'a = GoTc (V'/s + ^7) > (2-47) 

and 

{u^W)=0. (2.48) 

Since T does not depend on e or A, we can achieve the solution of Eq. (2.42) with 
A — ^ 00, by solving Eq. (2.47) without any limiting procedure. On the other hand, 
we have proven that Eq. (2.47) is equivalent to the Horiuchi's OCM Eq. (2.25), 
using Vb- Wc can thus prove the equivalence between Eq. (2.42) with A — >■ 00 
and Eq. (2.25); namely, the equivalence between the Kukulin's OCM and the 
Horiuchi's OCM, when the Pauli- forbidden state \u) is the exact eigen-state for 
the local potential Vb- In fact, this equivalence is also valid even when the \u) is 
not the eigen-state of Vb but the h.o. Pauli-forbidden state of the normalization 
kernel K, which is proved in the next subsection. 

For solving the 3a Faddcev equation, it is important to note that there exist 
some trivial solutions related to the orthogonality condition Eq. (2.21) for the T 
T-matrix. As is discussed in Ref. [9] in detail, the eigen-value solutions of the re- 
arrangement matrix S among the three different types of the Jacobi coordinates 
with the eigcn- value r = — 1 arc the redundant solutions of the Faddcev equa- 
tion Eqs. (2.17) or (2.26). For the h.o. Pauli-forbidden state \u), this eigen-value 
equation reads 

{u\S\ur) = r|D with r = -1 . (2.49) 
We find that the Faddeev component 



ro = Go\un 



with r = — 1 



(2.50) 



is a trivial solution with the total wave function I^'q = (1 + S)^jJQ = for r = — 1. 
The non-zero i/jq is the [21]-symmetric function with respect to the permutation 
of the 3a particles. We have two such trivial solutions for the 3a system with the 
total angular momentum L = 0. For this reason, the Faddeev equation (2.17) or 
(2.26) should be modified to^ 



AV' 



GoTS- ^ Go\un 



1 



{uriGown 



(2.51) 



in order to find a unique solution with A = 1. The solution of Eq. (2.51) with 
A = 1 automatically satisfies 



^l; = GoTS^, (u|(l + 5)1^^) = , {un 



for T 



-1 . 



(2.52) 



'^The inverse of the matrix elements in the last terms of Eq. (2.51) and Eq. (2.56) should be 
understood as the matrix inverse. 
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When the bound-state solution \ub) is used for \u), these r = — 1 eigen- values 
of Eq. (2.49) are no longer exactly r = — 1, but move to the r > — 1 values. The 
[3]-symmetric basis states constructed from 

^0= , / ={l + S)\uBn with T~-l (2.53) 
V3(l + r) 

involve some of the important shell-model like components such as |[3](04)) 
etc. [19] If these configurations are excluded, one cannot describe a compact shell- 
model like structure of ^^C. The original 3q: OCM equation should, therefore, be 
formulated by using 

P = \^o){%\ + P , (2.54) 

instead of P in Eq. (2.25). Since I^Z^o) involves a small admixture of the redundant 
components as 



{ubWo) = \n with r~-l , (2.55) 

the ground-state solution of the 3a system with the dominant components 
naturally involves a small admixture of the redundant components. We can also 
formulate an equivalent Faddeev equation to this modified 3q; OCM equation 
with P, which has a slightly different form from Eq. (2.51): 



ip . (2.56) 



The derivation of this equation and a detailed discussion of the almost redundant 
components of the Faddeev equation will be given elsewhere. [19] 

2.3 Equivalence between pairwise orthogonality conditions and the method of 
orthogonalizing pseudo-potentials in the 3-cluster systems 

In this subsection, we will prove the equivalence between the Horiuchi's OCM 
and the Kukulin's OCM, even when the Pauli- forbidden state \u) is not the eigen- 
state of the pairwise potential Vbj but the original h.o. Pauli- forbidden state of 
K. The essential point is that thejOCM T-matrix T{uj) defined through T{u!,e) 
in Eq. (2.24) is nothing but the T-matrix generated from the pseudo-potential 
in Eq. (2.30) with \ub) \u) and A ^ oo. 

Let us assume that \u) is not the eigen-state of Vb, but the h.o. eigen-state 
of K with the eigen- value 7 = 1. The T-matrix generated from V^^ with A — ^ 00 
is, from Eq. (2.31), 

f{u) = Tniu;) - {u; - Ho)G^{u;)\u) \ {u\Gd{u;){u; - Ho) , (2.57) 

where Td{lo) and Gd{oj) are the T-matrix and the full Green function of Vb, 
respectively. On the other hand, we separate V'^^^{e) in Eq. (2.23) as 

V^^^ie) = (e - Ho) - Aie -Ho- Vb)A = Vi + Xi\u){u\ (2.58) 
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with 

^1 = - {Ho + Vb)\u){u\ - \u){u\{Ho + Vd) , 

Xi = e+{u\Ho + VB\u) . (2.59) 

Then, the general formula in Appendix gives the T-matrix as 

T{u,e) = ti + (1 + hGo)\u) 1^^^^,^ I , H(1 + Goh) 

Ai - {u\Go + GotiGo\u} 

= ti + {u- Ho)Gi\u) ] A u\Gi{lo - Ho) , (2.60) 

Ai - {u\Gi\u) 

where ti and Gi are the T-matrix and the full Green function of Vi. Here we 
take e — oo (Ai — ^ oo) and obtain 

r(a;) = lim T{u;,e) = h - {u - Ho)Gi\u) . ,^ , A u\Gi{uj - Ho) . (2.61) 
£-►00 \u\G\\u) 

We note that uj — Hq — V\ \s expressed as 

u:-Ho-V^=uj-Ho-V^ + {Ho + Vj,)\u){u\ + \u){u\{Ho + Fd) 
= u- A{Ho + Vd)^ + \u){u\{Ho + V-d)\u){u\ 

= A{lo -Ho- Vb)A +\u)[cj + {u\{Ho + Vd)!^) ] {u\ . (2.62) 
Then, we find that Gi = (w — Ho — Vi)^^ is expressed as 

Gi = Ga + \u) J {u\ , (2.63) 

(jj + {u\Ho + Vt)\u) 

where Ga is defined by 

Ga = Gb- Gb\u) . , A u\Gji , (2.64) 

with Gd = {uj — Ho — Vr,)~^ , and satisfies 

A{uj -Ho- Vb)A Ga = Ga A{oj -Ho- Vb)A = A . (2.65) 
Since we have 

Gi\u) = \u) TTTT m— T etc. , (2.66) 

' ' ' ' Lo + {u\Ho + Vr>\u) ^ ' 

Eq. (2.61) becomes 

T(.) ff„)|„)_^^^J_^(„|(„ - ff„) . (2.67) 

On the other hand, Gi = Gq + Got\Go yields 

ti = (a; - Ho)Gx(}^ - Ho) - (u; - i?o) = - Ho)Ga{uj - Ho) -{to- Ho) 
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where the last term in the right-hand side cancels with the last term of Eq. (2.67). 
Thus, we find 

T{lo) = {lo- Ho)Ga{uj - Ho) - - Hq) 

= {u- Ho)Gb{u - Ho) -{u- Ho) -{u- Ho)Gb\u) . , . {u\Gb{uj - Ho) 
= Td - (o; - Ho)Gb\u) . . A u\Gb{u; - Ho) , (2.69) 

which is nothing but Eq. (2.57). We can also prove that the bound state pole 
of Vd is eliminated from Eq. (2.69), by using the spectral decomposition of Gd 
in Eq. (2.40). This proves the equivalence between the Kukulin's OCM and 
our Faddeev OCM, hence the equivalence between the Kukulin's OCM and the 
Horiuchi's OCM for the h.o. Pauli-forbidden state \u). 



3 3a OCM for the BFW aa potential 

As an example, we consider the 3a system interacting via the deep aa potential 
yBFW(^) = -122.6225 e-°-22 + 46^ erf (0.75 r)/r (MeV) (3.1) 

with r, being the relative coordinate between 2a's in fm. In Eq. (3.1), erf(a;) = 
(2/0F) /J^ e-* is the error function. The Pauli principle between 2a parti- 
cles are taken into account in terms of the bound states of F^^^(r). Since 
the Coulomb force is not exactly treated in the Faddeev formalism, we re- 
place the folded Coulomb term of Eq. (3.1) with the cut-off Coulomb force 



Table 1. Two-q (2a) bound-state energies of the BFW potential with and without the Coulomb 
force. The heading "momen. rep." stands for the 2a bound-state energies calculated in the 3a 
Faddeev code, by using the Born kernel including the cut-off Coulomb force with Rc = 10 fm. 
See the text for the momentum discretization points. RKG is the results by the Runge-Kutta- 
Gill method in the coordinate representation (full Coulomb). The Coulomb contribution is also 
shown in the rows "Coulomb cont." 



momen. rep. 


Os 


Is 


Od 


with Coulomb 


-72.6257 


-25.6186 


-22.0005 


without Coulomb 


-76.7051 


-28.7325 


-25.0908 


Coulomb cont. 


4.0794 


3.1139 


3.0903 


RKG 


Os 


Is 


Od 


with Coulomb 


-72.6255 


-25.6186 


-22.0005 


without Coulomb 


-76.7050 


-28.7325 


-25.0908 


Coulomb cont. 


4.0795 


3.1139 


3.0903 
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v^{r) = {Ae^ /r) 9{Rc — r) with Rc = 10 fm, introduced at the nucleon level. 
Here 9{x) is the Heaviside step function. For this type of the Coulomb force 
between protons, the aa Coulomb potential in Eq. (3.1) is replaced with 

yC^(r) = { erf (0.75 r) 

-(1/2) [erf (0.75 (r + Rc)) + erf (0.75 (r - Rc)) \] /r . (3.2) 

In the following, we use h^/Ma = 10.4465 MeV • fm^ and = 1.44 MeV • fm 
for the comparison with Ref. [15], unless otherwise specified. Table 1 shows the 
bound-state energies obtained by diagonalizing the negative-energy T-matrix 
in the momentum representation. We have two bound states, (Os) and (Is), 
in the relative S'-statc and one bound state, (Od), in the Z?-state, when the 
Coulomb force is included. The relative momentum p is descritized using the 
Gauss-Legendre ni-point quadrature formula for each of the four intervals of 
- 1 fm~^, 1-3 fm~^, 3-6 fm~^, and 6-12 fm~^. The large value of p beyond 
Pq = 12 fm~^ is also taken into account by using the Gauss-Legendre ria-point 
quadrature formula through the mappings = pQ-\-ia,n{'n{l + x)/4}.^ We choose 
ni = 15 and ns = 5, so that 65 points are used for solving the bound-state 
wave functions and the necessary T-matrix for solving the Faddeev equation. 
The second Jacobi coordinate q is also descritized with the Gauss-Legendre n2- 
point quadrature formula with n2 = 10 in the similar way to p, but in this 
case choosing only three major dividing points with q = 1 fm~^,3 fm~^ and 
6 fm^^ is good enough. The Gauss-Legendre ns-point quadrature formula with 
Q = + tan{7r(l -|-x)/4} is also applied to g > go = 6 fm~^. In Tabic 1, the 
result of the coordinate-space calculation using the Runge-Kutta-Gill method 
is also shown for comparison. In this case, wc use the full Coulomb force in 
Eq. (3.1). We find that the choice Rc = 10 fm in the momentum-representation 
calculation is accurate enough with the error less than 1 keV. This is also true 
even in the 3a Faddeev calculations. The aa phase shifts by the BFW potential 
in Eq. (3.1) is very well reproduced, as is shown in the original paper [14]. 

Table 2 shows the solutions of the 3a Faddeev equation (2.56) for the L = 
ground and excited states, obtained by using the BFW potential and the cut- 
off Coulomb force with Rc = 10 fm. The demarcation "b.s." implies that the 
bound-state wave functions of the BFW potential are used for the Pauli-forbidden 
states. Partial waves up to iuiax are included in 2a and {2a)-a relative motion. 
The convergence of the ground state energy is very rapid when we increase ^max 
from 4 to 10. The inaccuracy in ^max = 8 is within 1 kcV. Table 2 also shows 
the results of the variational calculations for Eq. (2.25) with Vd — > V^^^ and P 
being replaced by P in Eq. (2.54). Here we use the translationally invariant h.o. 
basis for the variational functions, [8, 9] with the maximum value of the total h.o. 
quanta A^max = 72. Agreement with the Faddeev calculation is satisfactory for 
the ground state. For the excited 0+ state, it deteriorates since A^max = 72 is still 
insufficient. We find a very large binding energy of 19.897 MeV for the 3a bound 
state, which is different from the result in Refs. [15] and [16]. This difference 
originates from the fact that we have used P instead of P. We have also calculated 

^These ns points for p are not included for solving the Faddeev equation (2.51) or (2.56), since 
these cause a numerical inaccuracy for the interpolation. 
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Table 2. Results of the 3a Faddeev calculations for the ground (Oi"'') and excited (O2"'') 0+ 
states of ^■^C, by using the BFW potential and the cut-off Coulomb force with Rc = 10 fm. The 
bound-state (b.s.) wave functions of the BFW potential, \ub), are used for the Pauli-forbidden 
states \u). Partial waves up to £max are included in 2a and {2a)-a relative motion. The heading 
ni-n2-n3 is the number of the momentum discretization points (see the text for details); "dim." 
stands for the full dimensionality of the diagonalizing matrix for the Faddeev equation (i.e., 
4ni(3n2 + n;j)(£max/2 + 1)); e-za is the expectation value of the 2a Hamiltonian with respect to 
the 3a bound-state solution; is the 3a bound-state energy; C(o4) is the overlap amplitude 
between the 3a bound-state wave function and the SU3 (04) configuration with the h.o. width 
parameter v — 0.28125 fm~^; and (/|/) in the last column is the squared norm of the redundant 
components. The results of the variational calculations, using the translationally invariant h.o. 
basis with the maximum h.o. quanta A'^max = 72, is also shown in the last rows for comparison. 



\u) 


^max 


711-712-^3 


dim. 


£2a (MeV) 


Esa (MeV) 


C(04) 


(/I/) 






4 


15-10-5 


6,300 


14.610 


-19.595 


0.9652 


2.7 X IQ- 


-4 


b.s. 


6 


15-10-5 


8,400 


14.485 


-19.894 


0.9613 


2.7 X IQ- 


-4 


O1+ 


8 


15-10-5 


10,500 


14.483 


-19.897 


0.9612 


2.7 X IQ- 


-4 




10 


15-10-5 


12,600 


14.483 


-19.897 


0.9612 


2.7 X IQ- 


-4 




h.o. 


variation 




14.482 


-19.897 




2.7 X IQ- 


-4 




4 


15-10-5 


6,300 


7.886 


-0.370 


0.1302 


1.9 X IQ- 


-6 


b.s. 


6 


15-10-5 


8,400 


8.431 


-0.485 


0.1410 


1.7 X IQ- 


-6 


O2+ 


8 


15-10-5 


10,500 


8.500 


-0.475 


0.1419 


1.7 X 10" 


-6 




10 


15-10-5 


12,600 


8.522 


-0.471 


0.1422 


1.7 X 10" 


-6 




h.o. 


variation 




9.950 


-0.241 




2.0 X 10" 


-6 



the overlap ampUtude of the ground-state wave function, C(o4) = {^l=o'^^\^)^ 
where the 5L'3 (04) sheh-model wave function is expressed in the 3a cluster 
model as [9] 

^L=o^^ = [^(40)(p)t^(40)(9)](o4)o = Y^R2o{p,bi)R2o{q,b2)Y(^oo)o{p,q) 

4 3 
-^;^^i2(p, bi)Ri2{q, 62)^(22)0 (P, q) + -^Ro4{p, bi)Ro4{q, b2)Y(^u)o{P, q) ■ 

(3.3) 

Here, ?7(4o)(p) and t/(4o)(q) are the single-particle SU3 states with the SU3 cou- 
pling (40) X (40) (04), Rne{x,v) is the radial part of the h.o. wave function, 
and ^(a£)l(P;9) the coupled angular-momentum function. The width parame- 
ters of the h.o. radial wave functions in the momentum representation are given 
by 61 = 1/47 and 62 = 3/I67 with = jiu = 2v. For the present calculation. 
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we have used the h.o. width parameter v = 0.28125 fm^^ in the configuration 
space.'^ Since the binding energy is very large, the C(o4) value is very close to 1, 
and is about 0.96. The squared norm of the redundant components defined by 
(/I/) with I/) = {ub\^) = {ub\{^ + are also shown in the last column in 

Tabic 2. 

Table 3 shows results of the 2>a Faddeev calculations when the Coulomb force 
is switched off. We find that the Coulomb contribution in the present 3a ground 

^This v value corresponds to a rather compact a cluster with the rms radius Va = (3/4y^) = 
1.414 fm. 



Table 3. The same as Table 2, but when the Coulomb force is switched off. In the upper half, 
denoted by "h.o.", the results for the h.o. Pauli-forbidden states \u) with the width parameter 
V = 0.28125 fm~^ are shown. 



\u) 




ni-n2-n3 


dim. 


e2a (MeV) 


E^a (MeV) 


C(04) 


(/I/) 






4 


15-10-5 


6,300 


12.405 


-25.296 


0.9740 


1.3 X 10" 


-12 


h.o. 


6 


15-10-5 


8,400 


12.262 


-25.563 


0.9707 


1.3 X 10" 


-12 


Oi+ 


8 


15-10-5 


10,500 


12.259 


-25.565 


0.9706 


1.3 X 10- 


-12 




10 


15-10-5 


12,600 


12.259 


-25.565 


0.9706 


1.3 X 10" 


-12 




h.o. 


variation 




12.258 


-25.565 


0.9706 








4 


15-10-5 


6,300 


6.836 


-6.181 


0.1343 


5.7 X 10" 


-12 


h.o. 


6 


15-10-5 


8,400 


7.117 


-6.404 


0.1423 


5.8 X 10" 


-12 


02+ 


8 


15-10-5 


10,500 


7.125 


-6.417 


0.1426 


5.8 X 10" 


-12 




10 


15-10-5 


12,600 


7.124 


-6.418 


0.1427 


5.8 X 10" 


-12 




h.o. 


variation 




7.165 


-6.412 


0.1430 








4 


15-10-5 


6,300 


12.722 


-27.431 


0.9738 


2.6 X 10" 


-4 


b.s. 


6 


15-10-5 


8,400 


12.582 


-27.745 


0.9700 


2.6 X 10- 


-4 


Oi+ 


8 


15-10-5 


10,500 


12.580 


-27.748 


0.9700 


2.6 X 10- 


-4 




10 


15-10-5 


12,600 


12.580 


-27.748 


0.9700 


2.6 X 10- 


-4 




h.o. 


variation 




12.579 


-27.748 




2.6 X 10- 


-4 




4 


15-10-5 


6,300 


9.062 


-5.731 


0.1121 


3.1 X 10- 


-6 


b.s. 


6 


15-10-5 


8,400 


9.327 


-6.060 


0.1213 


2.5 X 10- 


-6 


O2+ 


8 


15-10-5 


10,500 


9.329 


-6.077 


0.1217 


2.5 X 10- 


-6 




10 


15-10-5 


12,600 


9.328 


-6.078 


0.1217 


2.5 X 10- 


-6 




h.o. 


variation 




9.363 


-6.074 




2.5 X 10- 


-6 
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Table 4. Result of the 3a OCM with pairwise orthogonahty conditions (Horiuchi-type 3a 

OCM) for the 3a ground state, using the translationally invariant h.o. basis. The BFW potential 
for aa is used with H^/Mn = 41.7860 MeV • fm^ and v = 0.28125 frn'^ The Coulomb force is 
switched off. The heading A^max is the maximum value of the total h.o quanta included in the 
calculation; E2a is the lowest 2a diagonalization energy; e2a is the expectation value of the 2a 
Hamiltonian with respect to the 3a ground-state solution; E^a is the 3a ground-state energy; 
and C(o4) is the overlap amplitude between the 3a ground-state wave function and the SU3 (04) 
configuration. The Faddeev result in Table 3 with the h.o. |w) is also shown in the bottom row 
for comparison. 



TV 

-"max 


E2a (MeV) 


£2a (MeV) 


^3a (MeV) 


C(04) 


4 


9.81970 








6 


6.23729 








8 


2.34421 


14.0067 


-22.6177 


1 


10 


0.85837 


14.3470 


-22.6745 


0.9993 


12 


-0.11812 


12.9759 


-25.2357 


0.9811 


14 


-0.67019 


12.7748 


-25.2768 


0.9797 


16 


-1.03181 


12.4909 


-25.5086 


0.9740 


18 


-1.26800 


12.3918 


-25.5244 


0.9729 


20 


-1.43122 


12.3215 


-25.5528 


0.9716 


30 


-1.77280 


12.2608 


-25.5649 


0.9707 


40 


-1.86319 


12.2584 


-25.5652 


0.9706 


50 


-1.89424 


12.2583 


-25.5652 


0.9706 


60 


-1.90664 


12.2583 


-25.5652 


0.9706 


Faddeev 




12.2594 


-25.5653 


0.9706 



state with the BFW potential is 7.85 MeV, which imphes that our 3a ground 
state is rather compact compared with the microscopic 2>a cluster model. In 
the latter case, the standard value is 5 ~ 6 MeV. When the Coulomb force is 
neglected, we find that the second 0"*" state appears around E^^t ~ —6 MeV. In 
Table 3, we also show in the upper half, denoted by "h.o.", the results when 
the h.o. Pauli- forbidden states are used for \u), instead of the bound-state wave 
functions \ ub) of the BFW potential. The h.o. width parameter v = 0.28125 fm~^ 
is again used for this calculation. We find that the 3a ground state is less bound, 
but the energy difference is only 2 MeV. In this case, the elimination of the Pauli- 
forbidden states of the 3a system is rather easy, if we use the translationally 
invariant h.o. basis in the varaiational calculation. In Ref. [8], we have examined 
the equivalence between such a variational calculation and the present Faddeev 
calculation using the 2a RGM kernel. Tabic 4 shows the results of the 3a OCM 
with pairwise orthogonality conditions (namely, Eq. (2.25) with yBFW 
and P constructed from the h.o. |tt)), using the translationally invariant h.o. 



Y. Fujiwara, M. Kohno and Y. Suzuki 



17 



Table 5. The same as Table 4, but for the 3a OCM, using the Kukuhn's method of orthog- 
onahzing pseudo-potentials (namely, Eq. (2.42)) for the h.o. Pauli-forbidden states |m). The A 
parameter, A = 10^ MeV, is used to eliminate the h.o. Pauli-forbidden components. The paren- 
thesized numbers indicate the values for the modification of the last two digits when A = 10® 
MeV is used. 



-^max 


E2a (McV) 


e2a (MeV) 


(MeV) 


4 


9.81970 (56) 


6.25 X 10^ (6.25 x 10^) 


1.87 X lO'^ (1.87 X 10^) 


6 


6.23729 (08) 


3.02 X 10'^ (3.02 X 10^) 


9.07 X 10^ (9.06 X lO'') 


8 


2.34421 (01) 


14.0067 (69) 


-22.6177 (88) 


10 


0.85836 (16) 


14.3470 (77) 


-22.6745 (56) 


12 


-0.11812 (32) 


12.9759 (68) 


-25.2357 (70) 


14 


-0.67019 (37) 


12.7748 (58) 


-25.2768 (80) 


16 


-1.03181 (98) 


12.4909 (19) 


-25.5086 (98) 


18 


-1.26800 (17) 


12.3918 (29) 


-25.5244 (56) 


20 


-1.43122 (38) 


12.3215 (26) 


-25.5528 (40) 


30 


-1.77280 (95) 


12.2608 (19) 


-25.5649 (61) 


40 


-1.86319 (33) 


12.2584 (95) 


-25.5652 (64) 


50 


-1.89424 (38) 


12.2583 (94) 


-25.5653 (65) 


60 


-1.90664 (78) 


12.2583 (94) 


-25.5653 (65) 



basis. The BFW potential for aa is used with I^/Mn = 41.7860 MeV • fm^ and 
the h.o. Pauh- forbidden states with u = 0.28125 fm^^. The Coulomb force is 
switched off for simplicity. Since the present 3a ground state is very compact, 
the convergence with respect to the increase of A^max is very fast, and A^^max = 40 
is almost sufficient to obtain the converged result. 

We have also examined the 3a OCM, using the Kukulin's method of or- 
thogonalizing pseudo-potentials (namely, Eq. (2.42)) for the h.o. Pauli-forbidden 
states l^x). The results are shown in Table 5 with respect to the cases when the 
A parameter in Eq. (2.30) is A = 10^ McV and A = 10^ MeV. For A = 10^ MeV, 
we find a complete agreement with the results in Table 4 for N^nax ^ 8. When 
A = 10^ MeV is used, the energies deviate from the values in A = 10^ MeV in the 
last two digits, as long as the model space is large enough. This is, of course, a 
rather expected result, since the h.o. basis can conveniently eliminate the Pauli- 
forbidden states with some particular h.o. quanta within the finite number of 
the basis states. 

A conclusion derived from the above discussion is that the deeply bound 
feature of the 3a system in the BFW potential does not change appreciably even 
if one uses the real bound states for the BFW potential as the Pauli-forbidden 
states \u), as long as the dominant shell-model components are preserved in the 
Pauli-allowed space by using P. This is because the deeply bound states have 
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Table 6. The same as Table 4, but for the 2a folding potential obtained from the Schmid- 

Wildermuth force with v = 0.257 fm~^. As for the spin-isospin dependence, an almost pure 
Serber force with Xa = 2.4 is used. For the reduced mass, the standard value /Mn = 
41.4711 MeV • fm^ for the microscopic 3q cluster model is used. 



-^max 


E-2a (MeV) 


£2a (MeV) 


Eza (MeV) 


C(04) 


4 


14.3525 








6 


9.1136 


— 


— 





8 


4.9458 


17.6002 


-5.8189 


1 


10 


3.0476 


15.2471 


-6.5774 


0.9873 


12 


1.8698 


13.1271 


-9.1708 


0.9528 


14 


1.1573 


11.8925 


-9.6655 


0.9293 


16 


0.6854 


11.0975 


-10.1209 


0.9085 


18 


0.3637 


10.5721 


-10.3130 


0.8946 


20 


0.1348 


10.2268 


-10.4427 


0.8845 


30 


-0.3903 


9.6254 


-10.6095 


0.8661 


40 


-0.5598 


9.5415 


-10.6253 


0.8634 


50 


-0.6313 


9.5286 


-10.6271 


0.8630 


60 


-0.6665 


9.5264 


-10.6274 


0.8629 


72 


-0.6882 


9.5259 


-10.6274 


0.8629 



very large overlaps with the h.o. Pauli-forbidden states. The strongly attractive 
feature of the BFW potential in the 2>a system is related to the short range 
nature of this potential, in comparison with the usual folding potentials in the 
OCM formalism. As an example of the usual 3a OCM, we show in Table 6 
the result of the aa folding potential using the Schmid-Wildermuth force [20] 
with V = 0.257 fm"^. In this case, wc use H'^/Mn = 41.4711 MeV • fm^ for the 
microscopic 3a model. We find that the 3a energy is —10.63 MeV and it is not 
overbound, since the Coulomb energy is about 5-6 MeV. 

The short range nature of the BFW potential can be seen, by directly com- 



Table 7. Decomposition of the 3q energy E^a into the kinetic-energy and potential-energy 
contributions. Here the h.o. \u) is used. 





e2a (MeV) 


Ezoc (MeV) 


{Ho) (MeV) 


{V) (MeV) 


BFW 


12.258 


-25.565 


124.680 


-150.245 


SW 


9.526 


-10.624 


78.405 


-89.030 
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Table 8. The correlation between the 3a ground-state energy Esa and the depth of the direct 
potential Vq. For BFW the bound-state \ub} is used. 





Vo (McV) 


b (fm) 


E3a (MeV) 


SW (i/ = 0.257 fm-2) 


-97.7 


2.26 


-10.62 


SW {u = 0.275 fm-2) 


-105.4 


2.21 


-14.66 


BFW (b.s. \ub)) 


-122.6 


2.13 


-27.75 



paring the nuclear part of the aa potentials: 

yBFW ^ _ 122.6225 e-°-22 , 

Fd(SW) = -97.7 e-^-^'^^ . (3.4) 
In order to calculate the folding potential Vb (SW) , we have used the formula 

for the effective two-nucleon interaction v{r) = vqw e~'^'^ with w = W + BP^^ — 
HPt—MP^Pt-. We use an almost pure Serber force with the Majorana parameter 
m = 0.505, which corresponds to = + AB - AH - 2M = 2.4 in Eq. (3.5). 
The 1/e ranges of these potentials are h = 2.13 fm for the BFW potential and h = 
2.26 fm for the folding potential. Since we have calculated £20 , wc can evaluate 
the contributions of the kinetic-energy and the potential-energy terms separately 
through the simple expression; [Hq) = 2(3e2a — E) and (F) = 3(£' — 2£2a)- (See, 
for example, Ref. [21].) Table 7 clearly shows that the large 3a energy of the 
BFW potential is the result of the large cancellation of the kinetic-energy and 
potential-energy contributions. In the original 3a OCM calculation using the 
Schmid-Wildermuth force, Horiuchi [12] has used the value v = 0.275 fm~^ and 
the pure Serber force, = 2.445. In this case, the direct potential becomes 

Vd(SW) = -105.4 e-°-2°^ , (3.6) 

with b = 2.21 fm and the converged 3a bound-state energy is —14.66 MeV.^ 
We find that the 3a energy is strongly correlated with the depth of the direct 
potential Vq, as shown in Table 8. If we extrapolate the BFW value from the 
above two results of the SW force, we find —23.68 MeV for the BFW potential, 
which is close to the calculated value —27.75 McV. It is natural that a deep 
potential gives stronger binding in the 3a system, since the effect of the potential 
term is by factor 1.5 larger than in the 2a system, as was pointed out by Horiuchi 
[12]. 



"In Ref. [12] it is reported as -14.68 MeV. 
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4 Summary 

In this study, we have developed the Faddeev formahsm for the three cluster sys- 
tems, which "exactly" takes into account the Pauli-forbidden states. Actually, 
the exact Pauli-forbidden states of three-cluster systems are defined through the 
eigen-value problem of the three-cluster normalization kernel with the eigen- 
value zero. However, the pairwise orthogonality conditions to the total wave 
functions developed in this paper arc known to give a good approximation for 
the exact Pauli-allowed space obtained by the diagonalization procedure of the 
normalization kernel. For example, in the 3a system composed of the simple 
(Os)^ harmonic-oscillator (h.o.) shell-model wave functions, this correspondence 
is completely verified by enumerating the SUs allowed states in the translation- 
ally invariant h.o. basis. [22] 

The main result of this paper is that it is this type of three-cluster orthogo- 
nality condition model (OCM) with pairwise orthogonality conditions that leads 
to the complete equivalence to the three-cluster Faddeev equation interacting 
only by pairwise interactions. The pairwise interaction can be two-cluster RGM 
kernels with the linear energy dependence, two-cluster folding potentials of the 
effective two-body force, or the deep phenomenological local potentials like the 
Buck, Friedrich and Whcatley potential (BFW potential) [14]. In order to formu- 
late the three-cluster Faddeev equation with explicit elimination of the pairwise 
Pauli-forbidden components, we only need to use the modified T-matrix, T{io, s), 
which eliminates the off-shell singularity related to the existence of the Pauli- 
forbidden states. This T-matrix is, in general, energy dependent (e-dependcnt) 
for the three-cluster Faddeev equation using two-cluster RGM kernels, and the 
energy dependence is self-consistently determined by calculating the expecta- 
tion value of the two-cluster Hamiltonian with respect to the resultant Faddeev 
solution. [8] On the other hand, in the three-cluster OCM, considered in the 
present paper, T{uj) is energy independent. The Pauli-forbidden state in this 
OCM could be the h.o. wave functions of the microscopic two-cluster system, or 
the real bound states of the phenomenological local potentials between two clus- 
ters. However, the 3a Pauli-allowed space should be carefully defined, in order 
not to exclude the dominant SU^ components for the realistic description of the 
^^C ground state. We find that the Kukulin's method of orthogonalizing pseudo- 
potentials [13] is completely equivalent to the three-cluster OCM with pairwise 
orthogonality conditions. The latter was first proposed by Horiuchi [11]. We 
have proven this equivalence through the equivalence between either model and 
the present three-cluster Faddeev formalism using T[u)). A nice feature of the 
present Faddeev formalism is that the T-matrix description of the two-cluster 
interaction allows us to take the limit X ^ oo analytically in the method of or- 
thogonalizing pseudo-potentials, and that the solution of the Faddeev equation 
automatically guarantees the pairwise orthogonality conditions of the total wave 
function, owing to the orthogonality property of the T{u) T-matrix. 

As an example, we have applied the present three-cluster Faddeev formalism 
to the 3a system interacting via the BFW aa potential [14]. The Pauli-forbidden 
states are assumed to be the real bound states of the potential. We have found 
that this potential yields the 3a ground-state energy —19.897 MeV, which is a 
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different result from Refs. [15] and [16], in which the orthogonahty to the 2a 
bound-state solutions \ub) is very strictly demanded. This feature of the large 
overbinding does not change even when we use the h.o. Pauli-forbidden states 
with a reasonable width parameter of the a-clusters. This feature is traced back 
to the deep and short-range nature of the BFW aa potential, as compared with 
the folding potentials usually used in the microscopic 2>a cluster model. The 
validity of the So; boson model described by the BFW aa potential is examined 
by calculating other physical observables, like the rms radius of ^^C, using the 
obtained So; ground-state wave function. 
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A simple formula for T-matrix 

In this Appendix, we will show a general formula of the T-matrix for the surn of two potentials. 
From a simple calculation, the T-matrix for V = Vi -I- 1^2 is expressed as T = Ti -|- T2 with 

Ti=ti+ tiGoT2 , T2=t2+ tiGoTi , (A.l) 
where ti and t2 are the T- matrices of Vi and V2, respectively: 

ti^Vi+ ViGoti , t2 = V2 + V2Got2 ■ (A.2) 
If we further assume V2 = A |m)(m|, we find 

Thus, the solution of T = V + VGoT for V = Vi+ X\u){u\ is given by 

T = n + il+ ^igo)H ^_,_^^,^;^^^^^^^|^^ H(l + Got,) . (A.4) 

If we further move to A ^ cx), we find 

T = - (1 + ^iGo)|.)^^^^^^l^^(n|(l + GoU) , (A.5) 

with =V, + ViGoti. This expression leads to our basic relationship 

{u\ [1 + GoT\ = and [1 + TGo] |m) = . (A.6) 

As an example, let us use this formula in the method of orthogonalizing pseudo-potentials 

for Frgm(£). We set 

V'i^Vr,gm(£) and ti^T{uj,e) . (A.7) 

If we use the relationship between T{LO,e) and T(w,e), given in Eqs. (2.8), (2.9), and (2.10) of 
Ref. [8], we easily obtain 

T = T{uj,e)-[l + T{LO,s)Go{Lu)]\u) ^ 



(w|Go(w) -I- Go{Lj)T{u3,e)Go{oj)\u) 

x{u\ [l + Goiio)T{uj,e)] 

= T{cj,e)-{io-Ho)\u)-^{u\{LO-Ho)^f{uj,e) . (A.8) 

LO — £ 

Namely, our RGM T-matrix T{uj,e) is obtained from Vrgm{£) by the method of orthogonalizing 
pseudo-potentials. 
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